Majorana states and longitudinal NMR absorption in 3 He-B film 
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The topological superfluid 3 He-B supports massless Dirac spectrum of surface bound states which 
can be described in terms of the self-conjugated Majorana field operators. We discuss here the 
possible signature of surface bound states in nuclear magnetic resonance absorption spectrum in a 
3 He-B film. It is shown that transitions between different branches of the surface states spectrum 
lead to the non-zero absorption signal in longitudinal NMR scheme when the frequency is larger 
than the Larmour one. 

PACS numbers: 



I. INTRODUCTION 

Recently much attention has been devoted to the in- 
vestigation of bound fcrmion states on surfaces and inter- 
faces of topological superfluid 3 Hc-B. The presence of sur- 
face states in 3 Hc-B can be observed through anomalous 
transverse sound attenuation^— , surface specific heat 
measurements 5 -. Gapless fermion states are supported 
by the non-zero value of the topological invariant in 3 He- 
and have two - dimensional relativistic massless Dirac 
spectrum^ - — . Such massless fermions can be described 
in terms of the Majorana self-conjugated field operators 
which have been intensively studied recently in a number 
of condensed matter systems^. 

It has been suggested that NMR technique can be em- 
ployed to study the spectrum of surface states^ 3 -. The 
NMR measurements were implemented recently on su- 
perfluid 3 He films demonstrating various frequency shifts 
associated with the order parameter dynamics both at 
high-tcmperature A phase as well as with two non- 
equivalent low temperature B phase states characterized 
by the different values of the Ising variabl e 14 ' 15 . 

In this paper we focus on the contribution of the sur- 
face bound states to the ac magnetic susceptibility of 
He 3 B film. We demonstrate that the transitions be- 
tween different branches of surface states spectrum result 
in dissipation manifested in non-zero imaginary part of 
the Xzz component of magnetic susceptibility, where z is 
the axis normal to the film surface. This effect should 
provide the contribution to the longitudinal NMR ab- 
sorption signal absent both in the normal state and bulk 
He B phase. Under the action of the z component of 
magnetic field which destroys the self-conjugating Majo- 
rana states the magnetic absorption is suppressed at the 
frequencies smaller than the Larmour one. 



II. SPECTRUM OF SURFACE STATES IN 
3 HE-B FILM 



At first we introduce the basic formalism for treating 
the spectrum of fermionic quasiparticles. We consider 
3 He-B film confined in a slab at z > and homogeneous 
in x and y directions. The 3 Hc-B surface mode is derived 




from the quasiclassical BdG Hamiltonian 

H = -ihV z T 3 d z + fiA - |H • <x, (1) 

where A = Aijqi&j, and q = k/fcp where f,; are Pauli 
matrices of Bogolyubov-Nambu spin, a a are Pauli ma- 
trices of 3 He nuclear spin, 7 is the gyromagnetic ratio of 
the 3 He atom. 

The order parameter in 3 He-B is 3 x 3 matrix A a i where 
the Greek and Latin indices correspond to the spin and 
orbital variables. The geometrical confinement induced 
by the walls destroys the isotropy of the B phase order 
parameter—. As a result For the distorted 3He-B it has 
the form 



At 



where Aj_ is the gap for quasiparticles propagating in 
direction along the normal to the wall and Ay is the gap 
for quasiparticles propagating in directions parallel to the 
wall. 

Upon the specular reflection from the surface the quasi- 
particle momentum projection q z and therefore the part 
of the order parameter changes the sign which leads to 
the formation of surface bound states. To find the spec- 
trum of this states we employ the usual procedure, con- 
sidering the two parts of hamiltonian H = Hq + Hi so 
that 

H = -ihV z f 3 d z + hF(z) (2) 

T 

Hi = A||fi (a x q x + a y q y ) — -H • <x, (3) 

where F(z) = Aj_(z)a z q z and q = p/pf The hamil- 
tonian @ has zero energy eigenvalues corresponding to 
the degenerate surface bound states. To get into account 
correction from the perturbation terms Hi we find the 
eigenfunctions as a superposition 



vb = ^2Xj(pj(z), 



(4) 



i=i 
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where Xk arc the arbitrary coefficients and the generic 
terms ipj (z) are the eigenfunctions of hamiltonian H cor- 
responding to the zero energy e = 



i Pj (z) = A- 1 / 2 a j /3 j vq> \-K(z)}, 



(5) 



where A = (tpi\ifi) - 
cient and 

K(z) 



(ip2\f2) is the normalizing coeffi- 



Aj_ (s)ds. 
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The Pauli and Nambu spinors f3j and ctj satisfy the re- 
lations b z fi\.i — T0i,2 and f 2 ai i2 = ±0:1.2- 

Following the standard method we substitute the solu- 
tion in the form ((4]) into the equation (Ho + Hi)if> = eifj 
multiply by if>*(z) from the left and integrate over z. 
Then for the spinor X = (Xi,X2) T we obtain the two- 
dimensional Dirac equation 



[C&p + a z M]X = eX, 



(6) 



with the 'light velocity' given by 
1 f°° 

C=- — / dzAu(z)eiq>[-2K(z)} ~ AJp F 
Ap F Jo 

and 'mass' determined by the Larmour frequency M = 
huj H /2. 

For the massless particles we can choose the eigenfunc- 
tions of Eq. (IH) satisfying the relation 



X e = ia y X^ e 



(7) 



As a result the quasiparticle field operators are self- 
conjugated analogously to the Majorana fermions in rel- 
ativistic quantum field theories. In the external magnetic 
field the particles described by Eq.© become massive^ 
so that the property ([7]) does not hold and therefore 
the quasiparticlcs arc no longer self-conjugated Majorana 
fermions. 

The Dirac equation ((6]) determines the equation for the 
energy levels in the following form^ 



ei,2 = ±V(Cp) 2 + (hLo H ) 2 /4. 



(8) 



The energy spectrum of surface bound states ([8]) is sensi- 
tive only to the z projection of the magnetic field. It can 
result in a large anisotropy of magnetic susceptibility- 
if the magnetic field is much smaller then the effective 
dipolc field. However the larger magnetic field will reori- 
ent the spin axes eliminating the magnetic anisotropy^. 
Note that deriving the spectrum §E§ we have neglected 
the finite thickness of the slab. The size effect due to the 
overlap of quantum states localized at the opposite sur- 
faces of the slab leads to the splitting of Dirac coneii even 
in zero magnetic field. For sufficiently strong magnetic 
field this modification can be neglected. 



III. MAGNETIC SUSCEPTIBILITY AND NMR 
ABSORPTION SPECTRUM 

Now we consider the contribution of surface bound 
states to the imaginary part of ac magnetic susceptibil- 
ity component \ zz which determines the power absorp- 
tion under the experimental conditions of the longitudi- 
nal scheme of magnetic resonance when both the total 
magnetic field is directed along the z axis. 

To find the magnetic susceptibility let us use a conven- 
tional Kubo formula: 



Xij 



where 



T 
4 



d 2 p 

(2-kK)' 



Tr{a t G( P+ )a 3 G(p^}, (9) 



fc=l,2 



ioj M - £fc(p) 



is a temperature Green function and p± = (uj n ±e;/2, p). 
Here w n = ir(2n + 1)T is a fcrmionic and e/ = 27riT is 
a photonic Matsubara frequencies. The normalized wave 
functions if) 1,2 are given by the superpositions @ and 
correspond to the energy branches £1.2 in Eq.®. We use 
the value of the matrix element | (ipi \a z |-02) | = Cpje\ and 
the formula for the sum over fermionic frequencies 



T 



= /o(ei) - Mfg) 

- £i][ifi„_ - £2] iei-£ 2 +£i 



where fl n ± = w„ ± e;/2 and /0(e) = tanh(e/2T). Then 



changing the photonic Matsubara frequency by e; 
from Eq.(|9])we get that 

1 2 C 2 r p3dp (s 2 -e 1 )[f (e 1 )-f (e 2 )} 



8irh 2 



£ 2 {(friu)2 - (£2 - £l y 



iHui 



(10) 



where £1,2 = £1,2 (p)- We use now the dispersion relation 
(|HJ) and obtain that when the frequency is larger than the 
Larmour frequency uj > ojh the susceptibility given by 
Eq. fllOp has a non-zero imaginary part 



Imx z 



(—) 

Kmc) 



2 hu: 



-fa 



4T 



1 



H 



(11) 



Note that the non-zero dissipation Im\ zz 7^ occurs 
only when the frequency is larger than the threshold 
value (jjjj which is similar to the threshold behavior of 
absorption rate in semiconductors where the absorption 
edge frequency is determined by the band gap energy. 
However in the vicinity of the threshold the frequency 
dependence of absorption rate Eq.lfTTj) is completely dif- 
ferent from that of the electromagnetic wave absorption 
in the physics of semiconductors^. 

The estimate of imaginary magnetic susceptibility 
of a unit film area yields 



Imxz 



X"*A7° I It 



where Xn ~ J 2 kp/ Ep is a normal state susceptibility and 
£ = HVf/Ah is a coherence length. 
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IV. EFFECT OF SURFACE ROUGHNESS 

In general the quasiparticle energy levels are broadened 
due to the statistical fluctuation of the film surface which 
affect the boundary conditions for the wave functions. 
Different models of surface roughness related to the sur- 
face effects in 3 He were developed including the diffusive 
surface laye r 18 ! 19 , randomly rippled wall (RRW), ran- 
domly oriented mirrors (ROM) models 1 ^ and random 
scattering matrix model^i. Under the conditions of diffu- 
sive scattering the acoustic impedance data demonstrate 
the presence of surface bound states in 3 He-B film 2 -^ al- 
though there is no evidence of the relativistic massless 
Dirac spectrum. On the other hand the surface condi- 
tions can be varied in the experiments by coating the 
surface of several layers of 4 He22. For increased specu- 
larity factor the new features on the temperature depen- 
dence of acoustic impedance were observed indicating the 
formation of 2D Dirac energy spectrum^. 

Here we employ the ROM model assuming that the 
surface consists of small randomly oriented specularly 
scattering facets^. This model is applicable to describe 
the fluctuations with the scale much larger than kp 1 of 
the 4 He coated surface of the film. Within the frame- 
work of ROM model the important characteristic is the 
angle a s which constitutes the local normal vector to 
the wall n s with the z axis. Let us use the new co- 
ordinate system rotated by the angle a s with respect 
to the axis defined by v = z x n s . Then we obtain in 
the new coordinate system the expression for the order 
parameter matrix Ay = AikRkj, where R is the corre- 
sponding rotation matrix. Let us assume without loss of 
generality that the rotation axis coincide with the y axis. 
Then in the Eq.([5]) for the hamiltonian Hq we obtain 
F = (Aj_ cos a s a z — A|| sma s a x )q z and the perturbation 
term is given by 



Hi = T X q x (A||5- a cosa s + Aj_<7 z sma s 



7 



H-er. (12) 



To proceed further with analytical calculations we as- 
sume that the order parameter does not depend on the 
space coordinates so that Aj_,Aii = const. Then we 
obtain easily the zero energy eigenvectors of the hamil- 
tonian Hq in the form of Eq.([5|) with 



Ay sina. 



Aj_ cos a 8 + A 
A|| sin as 



1 



cos a s 



(13) 
(14) 



where A = (A|| sina s ) 2 + (Aj_ cosa s ) 2 . Correspond- 
ingly in the Eq.([5]) for the zero order wave functions we 
obtain 



K{z) 



1 



hv F 



Ads. 



The quasiparticle spectrum obtained along the perturba- 
tion theory scheme described above yields the spectrum 



in the form (|8]) but with modified parameters. We will 
study the modification of the absorption threshold which 
is determined by the 'mass' term and does not depend on 
the 'light velocity'. We therefore will neglect the modi- 
fication of 'light velocity' and focus on the 'mass' term 
which is given by 

(Aj_ cos a s + A) 2 - (Ay sino; s ) 2 
^- w ff( A±cosas+ A)2 + (A||sina 5 ) 2 ' ( ' 

To proceed further and calculate statistical average over 
the surface roughness we assume that \a s \ -C 1 and ob- 
tain to the leading order 



&H = 0J H 1 - 




(16) 



The above Eq. yields the fluctuating correction to the 
absorption edge in Eq. (|TT|) . It leads to the smoothing out 
the sharp absorption edge at the Larmour frequency. To 
estimate this effect we assume the Gaussian distribution 
of angle a s with the zero average value (a s ) = and the 
dispersion (a 2 ) = cr 2 . After that the average value of 
the susceptibility in the vicinity of Larmour frequency is 
given by 



where ao = 2 (Aj_/Ai 



) 2 1 — (lo/luh) 2 and 



S(ao, tr Q ) 



' exp [-{a/a a f 

(a - a )- 



-da (17) 



where a* = y^ao if «o > and a* = if ao < 0. 

The plots of the function (fPT|) for the different val- 
ues of a a are shown in Fig.([T]) demonstrating smoothing 
of absorption edge with increasing dispersion of surface 
ripples. Although in general for a a > the absorption 
signal is non-zero at the whole frequency domain however 
it is exponentially decaying for tu <C luh- The size of the 
crossover domain in Fig.([T]) is determined by the disper- 
sion 5ui = iono a . Therefore in general we can conclude 
that the absorption edge should well observed provided 
these fluctuations are small so that a a <C 1. 



V. CONCLUSION 

To conclude we have calculated the contribution of 
fcrmionic surface bound states to the ac magnetic sus- 
ceptibility of 3 He-B film. We have shown that in the 
longitudinal NMR scheme the non-zero absorption signal 
appears provided the frequency is larger than the thresh- 
old one determined by the Larmour frequency uj > lor- 
Such absorption is absent in the normal state of 3 He and 
can not occur due to the dynamics of the order parameter 
spin either. In zero magnetic field there is no frequency 
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FIG. 1: Plot of the function S(ao,a a ) for a a = 0.01 (solid 
line) and a a = 0.7 (dashed line). 



threshold for the dissipation which can be considered as 
the fingerprint of the gapless Majorana surface bound 
states. The surface fluctuations arc shown to smooth the 
threshold behavior out providing the small absorption in 
the frequency domain lo < luh- 
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